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3 '■ ABSTRACT 

>■ I We imaged the symbiotic star CH Cyg and two PSF cahbration stars using 

^^ I the unique 6.5 m MMT deformable secondary adaptive optics system. Our high- 

2 ; resolution (FWHM = 0.3"), very high Strehl (98 ± 2%) mid-infrared (9.8 and 

Q I 11.7 /xm) images of CH Cyg allow us to probe finer length scales than ever before 

^ I for this object. CH Cyg is significantly extended compared to our unresolved 

"^ ' PSF calibration stars (/i UMa and a Her) at 9.8 and 11.7 /im. We estimated the 

Qh' size of the extension by convolving a number of simple Gaussian models with 

O ' the /i UMa PSF and determining which model provided the best fit to the data. 

^ • Adopting the Hipparcos distance for this object of 270 pc, we found a nearly 

p3 '■ Gaussian extension with a FWHM at 9.8 /im of ~40.5±2.7AU (0.15±0.01") and 

.> ; a FWHM at 11.7/im of 45.9±2.7AU (0.17±0.01")- After subtracting out the 

/\ • Gaussian component of the emission (convolved with our PSF), we found a faint 

c^ • ~0.7" asymmetric extension which peaks in flux ~0.5" north of the stars. This 
extension is roughly coincident with the northern knotlike feature seen in HST 
WFPC2 images obtained in 1999. 

Subject headings: instrumentation: adaptive optics — binaries: general — stars: 
evolution — stars: formation — stars: Symbiotic — stars: individual (CH Cygni) 
— infrared: stars 



°The results presented here made use of the of MMT Observatory, a faciUty jointly operated by the 
University of Arizona and the Smithsonian Institution. 
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Introduction 



Symbiotic stars are a class of interacting binaries consisting of a hot component (usually 
a white dwarf) which ionizes the stellar wind of a cool component (generally a red giant or 
Mira, see Kenyon 2000). With typical component separations of ~1-50AU, these complex 
systems are responsible for some of the most dramatic Galactic jet sources (Eyres et al. 
2002). 

CH Cyg is one of the most enigmatic symbiotic stars - previous to 1963, CH Cyg showed 
no variable behavior and was actually used as a M6III spectral type calibrator (see Kenyon 
2001, for background details on CH Cyg's properties). Since 1963, CH Cyg has gone through 
multiple outbursting phases, producing bipolar radio jets in 1984 (Taylor et al. 1986), 1992 
(Karovska et al. 1998), and 1998 (Crocker et al. 2001), as well as X-ray jets (Galloway & 
Sokoloski 2004), and large ionized outflows (~10" in size, see Eyres et al. 2002; Corradi et 
al. 2001). B band flickering studies also point to the existence of an unstable accretion disk 
around the hot white dwarf component (Sokoloski & Kenyon 2003a,b). 

The CH Cyg system contains at least one red giant/white dwarf pair (Mikolajewski et 
al. 1987). Currently, there is ongoing debate (Sokoloski & Kenyon 2003b) as to whether the 
CH Cyg system is binary or the only known symbiotic triple star (consisting of either a red 
giant/white dwarf pair with a G star companion (Hinkle et al. 1993) or a red giant/white 
dwarf pair with an additional red giant companion (Skopal et al. 1996)). For the purposes 
of this paper, we adopt the interpretation of the CH Cyg system as an inner red giant /white 
dwarf pair with a period of 756 days (Hinkle et al. 1993) and an outer red giant with a period 
of 14.5 years (Skopal et al. 1996). All three stars are viewed nearly edge on and the outer 
giant periodically eclipses the inner binary. 

As an unresolved triple system, CH Cyg is very difficult to study. In order to completely 
understand the system it is necessary to disentangle the contribution of each component to 
the behavior of the system as a whole. Most studies of CH Cyg focus on the contribution of 
the hot white dwarf component to the system - jets (Galloway & Sokoloski 2004), ionized 
outflows (Eyres et al. 2002; Corradi et al. 2001), and hot component disk flickering (Sokoloski 
& Kenyon 2003a,b) primarily trace the properties of the hot component and the matter 
accreting on to it. In fact, the hot component dominates the total source properties at many 
wavelengths - to study the cool components, it is important to flnd a wavelength at which 
the cool components dominate over the hot component. To this end, we chose to study CH 
Cyg at mid-infrared (mid-IR) wavelengths, where the contribution of the hot component 
is negligible (<5% of the total emission, even during active periods). While near-IR and 
mid-IR photometry has previously been undertaken for CH Cyg (Taranova & Shenavrin 
2000, 2004; Bogdanov & Taranova 2001), CH Cyg has not previously been resolved at these 
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wavelengths, since sufficient spatial resolution at mid-IR wavelengths has not been available 
until recently (Close et al. 2003, Biller et al. 2005). Using the novel technique of adaptive 
optics at mid-IR wavelengths, we can directly image the emission from the cool component 
of the system by itself for the ffist time and study any circumstellar material on finer spatial 
scales than before. 



2. Observations and Data Reduction 

Data were taken on the night of 2003 May 13 (UT) at the 6.5 m MMT using the adaptive 
secondary mirror AO system (Wildi et al. 2003) with the BLINC-MIRAC3 camera (Hoffman 
et al. 1998; Hinz et al. 2000). Currently all other existing AO systems use a warm reimaged 
adaptive mirror and require several additional warm optical surfaces between the secondary 
and the detector. With the unique MMT adaptive secondary, we remove many of these 
warm optical surfaces from the optical path compared to a standard AO system. This 
results in a considerably reduced thermal background at the mid-IR, thus making adaptive 
optics possible for mid-IR observations. The adaptive secondary corrected the ffist 52 system 
modes and achieved Strehl ratios as high as 0.97±0.03 from 8.8 to 18 fira. These Strehl ratios 
are the highest ever presented in the literature for a large ground based telescope (Close et al. 
2003, Biller et al. 2005). For details on the MMT Adaptive Secondary AO system, see Wildi 
et al. (2003). We utilized the 128x128 SiAs BIB 2-20 /xm MIRAC3 camera (Hoffmann et al. 
1998). The 0.088"/pixel scale was used with the 9.8 and 11.7//.m 10% bandwidth filters. To 
remove thermal and detector instabilities we chopped at 1 Hz with an internal pupil-plane 
cold chopper in the interface dewar BLINC (Hinz et al. 2000) between the AO system and 
MIRAC3. 

We observed with the AO system locked continuously on CH Cyg. The 15° tilted BLINC 
dewar window is a high quality dichroic which refiected the visible light (A < 1 ^um) to the 
AO wavefront sensor and transmitted the IR through BLINC to MIRAC3. Since the internal 
chopper in BLINC was past the dichroic, continuous 1 Hz chopping did not affect the visible 
light beam and hence the AO lock was unaffected. To further calibrate the background (in 
addition to chopping) we nodded ~ 6-8" in the telescope's azimuth direction (the horizontal 
direction in Figure 1) every minute. The internal chopper was set to run in the altitude 
direction (the vertical direction) with a chop throw of ~20". The derotator was disabled 
during these observations to help minimize the residual background structure. 

The 0.505" PA = 269° Washington Double Star catalog astrometric binary WDS 02589+2137 
BU was observed during a later run (2003, November 25 UT) and used to calibrate the cam- 
era's orientation and its 0.088"/pixel platescale. 
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For the 9.8 and 11.7 ;um filters we obtained 2x1 minute coadded chop differenced images 
(one image from each nod) for CH Cyg. We utihzed a custom IRAF script to reduce this 
mid-IR data (Biller et al. 2005). The script produced two background subtracted images by 
subtracting the image taken at the first nod position with that taken at the second, and vice 
versa. Images were then bad pixel corrected, flat-fielded, cross-correlated, and aligned to a 
reference nod image (to an accuracy of ~0.02 pixels). Mid-IR images of the PSF calibration 
stars (/i UMa and a Her, both of these stars are binary but appear unresolved to MIRAC, 
see Close et al. [2003] for a discussion of PSF suitability) were obtained and reduced in an 
identical manner to CH Cyg. Reduced images of CH Cyg and the two PSF stars are presented 
in Figure 1. In order to directly compare our science images with the PSF images, we have 
not rotated the images to place north up - azimuth is vertical and altitude is horizontal. At 
the time of the observations, the Parallactic Angle (PA) was roughly 160°. 

The CH Cyg images appear slightly yet significantly extended compared to images of 
the two PSF stars. To highlight this extension, we present one-dimensional normalized cut 
plots in both altitude and azimuth directions across the peak for the 9.8 and 11.7/im images 
of CH Cyg and PSF star \x UMa in Figures 2 and 3. Three rows/columns were averaged 
together for each cut plot to increase signal to noise. Since CH Cyg is a very bright object 
(565 Jy at 12 /xm) and errors are dominated by photon noise, error bars for these plots are 
considerably smaller than the size of the plot symbols and have not been plotted. In all 
the altitude and azimuth cutplots at 9.8 and 11.7/im, the CH Cyg profiles are noticeably 
extended compared to those for the /x UMa PSF. We measured the FWHM, ellipticity, and 
positional angle of any such ellipticity for all the images in Figure 1, as well as for images 
of RV Boo, an AGB star with a small extended dust disk observed the same night (Biller 
et al. 2005). In Figure 4 we plot the ellipticity and FWHM for our datasets. CH Cyg has 
a similar eccentricity but a significantly higher FWHM than the PSF stars. Thus, CH Cyg 
looks clearly extended compared to the PSF stars. 

In order to detect the gas counterpart to the mid-IR emitting dust, CH Cyg was observed 
in the CO (1-0) transition at the Arizona Radio Observatory 12 m telescope on Kitt Peak 
on the night of 29 June 2005. No CO (1-0) was detected in 4 hours of integration time. 



3. Analysis 

In order to study small-scale circumstellar structure around CH Cyg, we deconvolved 
our high S/N (S/N > 1000) images of CH Cyg using both the PSF stars /x UMa and a Her 
(see Biller et al. 2005, for details on PSF star selection) with the IRAF Lucy deconvolution 
task (Lucy 1974). So as not to introduce noise amplification artifacts due to deconvolution. 
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we deconvolved images for only 10 iterations. To test whether artifacts were being introduced 
by the Lucy algorithm, we also deconvolved images of the /x UMa PSF with those of the 
a Her PSF in order to establish the expected result of deconvolving one point source by 
another point source. Deconvolved images of CH Cyg and the PSF stars appear in Figure 5. 
Due to the very high Strehl (~0.98) and high signal-to-noise ratio in our PSF images we 
could detect low contrast structure on "super-resolution" scales of ~0.2" with the IRAF 
Lucy deconvolution task (Biller et al. 2005). After deconvolution, CH Cyg still appears 
significantly extended compared to the deconvolved PSF images. 

To estimate the size of the extension around CH Cyg, we convolved a range of simple 
Gaussian dust envelope models (with FWHMs of 0.1, 0.14, 0.15, 0.16, 0.18, and 0.2" before 
convolution) with the /i UMa PSF. We determined the best fit model by comparing model 
cut plots with data cut plots; the best fit model minimizes the standard deviation of [cut- 
plotiddXd)- cutplotimodel)]. At 9.8 /im, the 0.15" model provided the best fit to the altitude 
cut plot, while the 0.14" model provided the best fit to the azimuth cut plot. At 11.7/im, the 
0.17" model provided the best fit to the altitude cut plot, while the 0.16" model provided the 
best fit to the azimuth cut plot. Best fit model cut plots are presented in Figures 2,3. The 
slightly broader cut plot profile along the altitude direction in both these cases is probably 
due to instrumental effects and is not physical. Utilizing the Hipparcos distance of ~ 270 pc, 
we adopt a FWHM at 9.8 /xm of -40.5±2.7AU (0.15±0.01"), and a FWHM at 11.7/xm of 
45.9±2.7AU (0.17±0.01"). 

We scaled and subtracted the best fit Gaussian model (convolved with the /i UMa PSF) 
from the data in order to reveal any asymmetric structure. Images of residual asymmetric 
structure are presented in Figure 6. The residual structure peaks in intensity ~0.5" north 
of the star (marked by a diamond) and is roughly coincident with the extension seen by 
Crocker et al. (2001) at ~0.7" north of the star in HST WFPC2 images (marked by the 
letter A in Fig. 6). The structure also appears to be extended along a similar axis to the 
radio/optical jets. 



4. Modeling 

Since the residual asymmetric structure we observe around CH Cyg is faint compared to 
the symmetric Gaussian structure (peak fiuxes only ~ 2% that of the symmetric component) 
we have concentrated our modeling efforts on the symmetric component only. We model the 
IR emission of CH Cyg as an optically thin shell passively heated by the star. 

We approximate the stellar photosphere by the Kurucz (1979) model spectrum for 
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M6IIIe stars with an effective temperature of Tefj = 3000K. The dust is taken to be amor- 
phous sihcate since CH Cyg is a M star. Taking all grains to be spherical with a being 
their spherical radius, we assume a power-law dust size distribution dn{a)/da oc a~" which 
is characterized by a lower-cutoff flmin, upper-cutoff Cmax and power-law index a. For sim- 
plicity, we take amin = 0.01/im since smaller grains will undergo single-photon heating, and 
Omax = 1000 /im since larger grains are not well constrained by the currently available IR 
photometry. 

We approximate the dust spatial distribution with a modified power-law dn/dr oc 
(1 — rmin/r) (rmin/'r)'' where rmin is the inner boundary of the envelope which we take to be 
the location where silicate dust sublimates. For CH Cyg with a luminosity L^ ^ 6900 Lq, 
submicron-sized silicate dust achieves an equilibrium temperature of T ^ 1500 K and starts 
to sublimate at r ^ 2AU. Therefore we take rmin = 2AU. This functional form has the 
advantage that on one hand, it behaves like a power-law dn/dr oc r^'^ at larger distances 
(r ^ rinin), and on the other hand it peaks at r^ = r^i^ {P + 7) /7, unlike the simple power- 
law which peaks at rmin. The latter is unphysical since one should not expect dust to pile 
up at rmin where dust sublimates! We take 7 = 2 as expected from a stationary outflow with 
a constant mass loss rate (e.g. see Habing, Tignon, & Tielens 1994). But other 7 values 
have also been considered (see below), in view of the time- variable nature of this object. We 
take Tmax = 200 AU which is large enough for our IR emission modeling purpose since there is 
very little dust beyond this outer boundary which emit at the wavelengths we are interested 
(A < 200 ;um). Therefore, we have only two free parameters: a - the power-law exponent for 
the dust size distribution and j3 - the dust spatial distribution parameter which determines 
where the dust peaks. 

Using the dielectric functions of "astronomical silicates" (Draine & Lee 1984) and Mie 
theory (Bohren & Huffman 1983), we calculate the absorption cross sections of amorphous 
silicate grains as a function of size, as well as their equilibrium temperatures as a function of 
radial distance from the central star (in thermal equilibrium with the illuminating starlight). 
We then obtain the dust model IR emission spectrum by integrating the dust emission over 
the entire size range and the entire envelope. 

The model spectrum is then compared with (1) the 9.8 /im and 11.7/im photometry from 
our MMT AO images, (2) the 12, 25, 60, and 100 /xm IRAS {Infrared Astronomical Satellite) 
broadband photometry, (3) the 7.7-22.7 /im mid-IR low-resolution spectroscopy obtained 
with the IRAS Low Resolution Spectrometer (LRS; with a resolution of A/ A A ~ 20), (4) the 
2.5-45 /im near- to mid-IR high-resolution spectroscopy obtained with the Short Wavelength 
Spectrometer (SWS; A/AA ~ 2000) on board the Infrared Space Observatory (ISO), and (5) 
the 43-197 /xm mid- to far-IR medium-resolution spectroscopy obtained with the ISO Long 
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Wavelength Spectrometer (LWS; A/AA ~ 150-200). Given that CH Cyg is a variable star, it 
is to our great surprise that these observational data taken at different time epochs are in 
reasonably good agreement with each other: the MMT-AO and IRAS fluxes are only ~ 12% 
and ~25%, respectively, smaller than the ISO SWS flux (see Fig. 7). The SWS spectrum 
is also in close agreement with the K, L, and M photometry. But the IRAS LRS flux is 
lower than the SWS spectrum by a factor of ~ 16. We therefore take the ISO SWS and LWS 
spectra as the comparison basis and increase the MMT-AO fluxes, the IRAS broadband 
fluxes, and the IRAS LRS spectrum by a factor of 1.12, 1.25 and 16, respectively. This 
brings the MMT-AO and IRAS data into agreement with the ISO data. 

The model which best fits the observed IR emission spectrum and photometric points 
has a Ri4.1 and (3 f^^lO (see Fig. 7), with the dust distribution peaking at r^ ^ 12 AU, a total 
visual optical depth Ty ~ 0.16, and a total dust mass rridnst ~ 1-04 x 10^'^ g ~ 5.2 x 10^'^ M0. 
It is not easy to translate the total dust mass into the mass loss rate which is actually 
considered as one of the most important tasks in the study of stellar evolution (e.g. see 
Willson 2000). To estimate the mass loss rate of CH Cyg, we need to assume a gas-to-dust 
mass ratio rrigas/'n^dust and a gas outflow velocity v. With mgas/'n^dust = 100 and w = 50 kms~^ 
(based on the velocity of the Ca-II absorption component which, according to Skopal et al. 
1996, is likely associated to the expanding circumstellar wind), we estimate the mass loss 
rate of CH Cyg to be ^ 2.7 x 10"^ M© yr'^. 

In order to compare the best fit model directly to our MMT AO images, we convolved 
the model radial profile with our AO PSF image of /i UMa. After convolution, the best 
fit model to the observed IR emission spectrum and photometric points produces an image 
appreciably larger than the MMT AO image. We have considered models with a grid of 
parameters (/3 = 2-20, 7 = 1-4, r^ax = 25-250 AU). While multiple models fit the IR emission 
spectrum reasonably well, none of the models is able to simultaneously replicate the AO 
image and the observed IR emission. For instance, the best fit model to the IR emission 
spectrum {a ?^4.1 and (3 f^ilO) has a FWHM ~ 0.4" at 9.8 jum (after convolution with the 
AO PSF), while the AO image of CH Cyg has a FWHM < 0.37". Given the extremely high 
Strehl ratios during these observations (~100%) and the resultant high image fidelity, this 
is a significant divergence! These models are likely too simple to fully describe this complex 
system. The models assume only one star, with a steady state outflow. However, the CH 
Cyg system contains 3 stars - we have not taken into account the disruptive effects upon the 
dust of the other two stellar components of this system. The assumed functional form for 
the dust spatial distribution may be too simplified. We should also note that, as a variable 
star, the spatial and size distributions of the dust (which determine the surface brightness 
distribution of the dust shell around CH Cyg) at the epoch when the AO images were taken 
can be quite different from those when the IRAS photometry and the ISO spectroscopy were 



obtained. As a result, it is very unlikely that a dust model will fit both the observed IR 
emission spectrum and the AO images. Thus this is not a test of the validity of the model. 



5. Discussion 

Circumstellar dust in the CH Cyg system has previously been studied by Taranova 
and collaborators. With near- and mid-IR photometric observations, Taranova & Shenavrin 
(2004) follow the properties of a circumstellar dust envelope around CH Cyg. According 
to Taranova & Shenavrin (2004) 's model, mass was injected into the cool component's cir- 
cumstellar dust envelope around 1983 by an as yet unknown mechanism, coincident with 
the bipolar radio outfiow observed by Taylor et al. (1986). The optical depth of the dust 
envelope increased through 1992. By 1996, the dust envelope began to dissipate and had 
decreased to pre-mass injection optical depths by 2001. 

Taranova & Shenavrin (2000, 2004) model the circumstellar dust envelope of CH Cyg 
as a single temperature shell around the cool component which varies in size and total dust 
mass during mass injection episodes. From 1983 to 2003, they claimed that the dust shell 
radius varied between ~ 10^^-10^^ cm, or from about 7-70 AU (but Bogdanov & Taranova 
[2001] argue for a much large extension: 15-15000 AU). We observe CH Cyg to be similarly 
extended, with a FWHM of ~ 40-46 AU. 

The symmetric extended dust structure around CH Cyg is considerably larger than the 
orbits of the stellar components. In order to estimate the semimajor axis of the inner 756 day 
period binary, we adopt a white dwarf/hot component mass of ~ 0.2 Mq and a red giant/cool 
component mass of ^2Mq (Hinkle et al. 1993). We calculate a semimajor axis estimate of 
~2AU. To estimate the semimajor axis of the outer 14.5 year period binary, we adopt an 
outer red giant mass of ~ 2 Mq and find a semimajor axis estimate of ~ 10 AU. The observed 
dust envelope extension of ~ 40-46 AU is larger than both the inner and outer binary orbits. 
Thus, the dust envelope contains all three stellar components of this system. The fact that 
all three stellar components of this system lie within the dust envelope may be significant 
in explaining why our AO dust shell images diverge from our spherically symmetric dust 
envelope models - possibly, the stellar components of this system have a hand in shaping the 
dust envelope into something other than a classic, red giant wind spherical shell. However, 
whatever specific contribution the stellar components have had on the shape of the dust 
envelope remains unresolved and unquantified. 

After subtracting out the symmetric component of the emission, the residual emission 
appears extended along a similar north-south axis to those of the radio jets reported by 
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Crocker et al. (2001). The residual structure peaks in flux ~0.5" north of the star and 
is roughly coincident to the knotlike structure observed by Crocker et al. (2001) at 0.7" 
north of the star in HST WFPC2 images (marked by the letter A in Fig. 6). The structure 
also appears to be extended along a similar axis to the radio/optical jets, and is similar to 
the extent of the radio jets and ionized optical structure. However, while bipolar structure 
is observed in both the radio and optical, the mid-IR structure is unipolar. One possible 
explanation is that mid-IR emitting dust is denser towards the northern lobe of the bipolar 
outflow. The outflow shocks and heats the dust towards the north, hence a bright mid-IR 
dust structure is observed. Towards the south, the outflow is unimpeded by dust and, thus, 
has a larger optical and radio extent, but no mid-IR counterpart. 



6. Conclusions 

The combination of adaptive optics with a deformable secondary produces very high 
Strehl images in the mid-IR and high PSF stability regardless of the seeing, airmass, or 
target brightness. Such PSF stability allows us to determine object properties with very 
high fidelity. Adopting the Hipparcos distance to CH Cyg of ~ 270 pc, we found a nearly 
spherically symmetric extension with a FWHM at 9.8 ;um of ~40.5±2.7AU (0.15±0.01") 
and a FWHM at 11.7/im of ~45.9±2.7AU (0.17±0.01"). After subtracting out a model 
of the spherically symmetric emission, we found residual extension ~0.5" north of the star 
roughly coincident with the optical knotlike feature observed by Crocker et al. (2001). 

Adaptive optics at mid-IR wavelengths is a very promising new technique that allows 
for uniquely stable PSFs and high Strehls. A high degree of PSF stability will eliminate 
morphological ambiguities due to poor (seeing-limited) PSF calibrations. Mid-IR AO should 
have a significant impact on any field where mid-IR imaging is possible. 

These MMT observations were made possible with the hard work of the entire Center for 
Astronomical Adaptive Optics (CAAO) staff at the University of Arizona. Graduate student 
Wilson Liu helped run the MIRAC3 camera during the run. The adaptive secondary mirror 
is a joint project of University of Arizona and the Italian National Institute of Astrophysics 
- Arcetri Observatory. We would also like thank the whole MMT staff for their excellent 
support and flexibility during our commissioning run at the telescope. 

The secondary mirror development could not have been possible without the sup- 
port of the Air Force Office of Scientific Research under grant AFOSR F49620-00-1-0294. 
BAB is supported by the NASA GSRP grant NNG04GN95H and NASA Origins grant 
NNG056L71G. LMC acknowledges support from NASA Origins grant NNG056L71G and 
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NSF CAREER award. AL acknowledges support from the University of Missouri Sum- 
mer Research Fellowship, the University of Missouri Research Board, and the NASA award 
P20436. JHB acknowledges support from NSF grants AST-9987408 and AST-0307687. 
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Fig. 1. — The 9.8 and 11.7 /xin images of CH Cyg and PSF stars fi UMa, a Her, and AC Her 
as observed at the MMT. The first, second, and third Airy rings are visible in these images. 
CH Cyg appears shghtly extended compared to the PSF stars - while the PSF stars show 
a peak distinct from the first Airy ring, the slightly extended emission from CH Cyg blurs 
into the first Airy ring. The faint point source in the far lower left of each MMT image is a 
MIRAC3 ghost. 
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Fig. 2. — One dimensional cut plots along the altitude (top) and azimuth (bottom) directions 
for the 9.8 yum images of CH Cyg and the PSF star /i UMa. The CH Cyg data is plotted as 
crosses, the /x UMa PSF is plotted as a dotted blue line, and the best fit Gaussian model 
convolved with the n UMa PSF is plotted as a solid red line. Error bars are small compared 
to the size of the plot symbols. CH Cyg appears slightly but notably extended relative to 
the /i UMa PSF. 
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Fig. 3. — One dimensional cut plots along the altitude (top) and azimuth (bottom) directions 
for the 11.7 /im images of CH Cyg and PSF star fi UMa. The CH Cyg data is plotted as 
crosses, the /x UMa PSF is plotted as a dotted blue line, and the best fit Gaussian model 
convolved with the n UMa PSF is plotted as a solid red line. Error bars are small compared 
to the size of the plot symbols. CH Cyg appears slightly but notably extended relative to 
the /i UMa PSF. 



-15- 



0.370 


* 


1 ■ ■ ■ 


; 






i 


^-^ 


CH Cyg 




z 


^ 0.360 


/ ij 




-_ 


w 


'. 




z 


o 


I 




z 


V. 


- 




- 


^ 0.350 


L 


RV Boo 
* 


_I 


E 


; 


z 


:l 


; 




': 


°°. 0.340 


- 




- 


CD 


: 




: 


@ 


va Her 




: 


E 0.330 


XX 




-; 


sz 


- 




- 


5 


XM UMo 




': 


0.320 


r ^'AC Her 




-. 




1 . . . 


1 . . . 


- 



0.00 



0.02 
eccentricity 



0.04 
9.8 /im 



0.06 



0.410 






* 








^ 0.400 

w 

o 

^ 0.390 


j 


Ch 


Cyg 






j 


E 

1^ 0.380 

* 


i- 










■\ 


@ 0.370 


:_ 






X 


X 


_: 


£ 


: 






AC Her 


a Her 


: 


"i 0.360 

M— 


~ 


M 


UMo 






~ 


0.350 















0.05 0.06 0.07 0.08 0.09 

eccentricity @ 11.7 /^m 



0.10 



Fig. 4. — The 9.8 yum (top) and ll.T/iin (bottom) FWHM and ellipticity of CH Cyg and the 
PSF stars /i UMa and a Her. CH Cyg appears shghtly extended in FWHM compared to 
the PSF stars. 
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Fig. 5. — Deconvolved images of CH Cyg. The plate scale (after magnification) is 
29.3 mas pixel"^. Deconvolved images of the n UMa PSF (deconvolved using a Her as a 
PSF) is shown on the right for comparison. CH Cyg appears slightly extended after decon- 
volution compared to the PSF stars. 
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Fig. 6. — Top: Faint residual structure (<2% of the peak) left over after scaling and sub- 
tracting the best-fit Gaussian model (convolved with the /j, UMa PSF) from the data. North 
is up and east is left in these images. The diamond point represents the location of the stellar 
components of the CH Cyg system in these images. Bottom: An HST WFPC2 image of CH 
Cyg from August 1999 (Crocker et al. 2001) is shown for comparison. The residual extension 
in the mid-IR is roughly coincident with the northern knotlike/bipolar feature observed by 
Crocker et al. (2001) in the optical and the radio. 
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Fig. 7. — IR emission and model fit to the CH Cyg dust shell. Open black diamonds - U, 
B, V, R, J, H, K, L, and M photometry from literature; open black triangles - MMT-AO 
photometry at 9.8/im and 11.7/im (increased by a factor of 1.12); open black squares - 
IRAS photometry at 12, 25, 60 and 100 /xm (increased by a factor of 1.25). Blue dashed line 
- stellar photospheric spectrum; solid blue line - ISO SWS spectrum; solid cyan line - ISO 
LWS spectrum; solid magenta line - IRAS LRS spectrum (increased by a factor of 16). The 
dust model is shown as a green dotted line, while the solid red line plots the sum of the dust 
and the star. 



